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Effect of the band structure topology on the minimal conductivity for bilayer graphene with 

symmetry breaking 
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Using the Kubo formula we develop a general and simple expression for the minimal conductivity in systems 
described by a two by two Hamiltonian. As an application we derive an analytical expression for the minimal 
conductivity tensor of bilayer graphene as a function of a complex parameter w related to recently proposed 
symmetry breaking mechanisms resulting from electron-electron interaction or strain applied to the sample. 
The number of Dirac points changes with varying parameter w, this directly affect the minimal conductivity. 
Our analytic expression is confirmed using an independent calculation based on Landauer approach and we find 
remarkably good agreement between the two methods. We demonstrate that the minimal conductivity is very 
sensitive to the change of the parameter w and the orientation of the electrodes with respect to the sample. Our 
results show that the minimal conductivity is closely related to the topology of the low energy band structure. 
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Introduction. — After the first quantum Hall measurement 
on graphene 111 \m the physics of graphene has become one 
of the leading research field in physics. The bilayer graphene 
has been studied first experimentally [3J by Novoselov et al. 
and theoretically Jj] by McCann and Fal'ko. Recent observa- 
tions [5, 6] indicate that in bilayer graphene spontaneous sym- 
metry breaking may arise from electron-electron Coulomb in- 
teractions. Occurance of such broken-symmetry states gener- 
ates more attention to the topological changes in Fermi surface 
in high quality suspended bilayer graphene devices. Lemonik 
et al. studied the spontaneous symmetry breaking and Lif- 
shitz transition in bilayer graphene [7]. Spontaneous inversion 
symmetry breaking in graphene bilayers has also been inves- 
tigated by Zhang et al. ||8|]. Vafek and Yang applied renor- 
malization group approach to study the many-body instabil- 
ity of Coulomb interacting bilayer graghene [9J. Nandkishore 
and Levitov JlOll . and Gorbar et al. jllll studied competition 
between different ordered states in bilayer graphene. Sponta- 
neous symmetry breaking in two-dimensional electronic sys- 
tems with a quadratic band crossing was studied by Sun et 
al. II12II . The quantum theory of a nematic Fermi fluid has 
been proposed by Oganesyan et al. J13J]. 

The low energy Fermi surface of bilayer graphene is dom- 
inated by the trigonal warping, first shown by McCann and 
Fal'ko |4]. At very low energies, namely below the Lifshitz 
energy the trigonal warping results in a breaking of the con- 
stant energy lines into four pockets. The Lifshitz energy is 
typically of order of 1 meV lil4]. Besides many other effects 
the trigonal warping has a substantial influence on the mini- 
mal conductivity cr™'". 

Without trigonal warping in Ref. Il5| it was found that 
cr™° = 8 (To, where ctq — e'^/inh). This later was confirmed 
by Snyman and Beenakker 1II6II using the Landauer approach. 
Actually, the calculated minimal conductivity can take non- 
universal values depending on the order of the dc limit and 
the integration over energies as shown by Ziegler JlTll . The 
importance of the order of the frequency and the temperature 



limit was pointed out by Ryu et al. OlSll . Trushin et al. showed 
that electron-hole puddle formation is not a necessary condi- 
tion for finite conductivity in bilayer graphene at zero aver- 
age carrier density 1.19.1 . Culcer and Winkler studied the role 
of the external gates and transport in biased bilayer graphene 
using the density operator formalism and quantum Liouville 
equation 12011 . From self-consistent Born approximation 
Koshino and Ando found that in the strong-disorder regime 
cr™" = 8 (To. while in the weak-disorder regime ct™" = 24 o-q- 
In Ref. [21.1 the role of the trigonal warping was studied and 
it was shown that the contributions of the four pockets to the 
minimal conductivity gives cr™'" = 24 ao- Moghaddam and 
Zareyan showed that the minimal conductivity of graphene 
bilayers is anisotropic with respect to the orientation of the 
connected electrodes when the trigonal warping is taken into 
account 1221 . 

The spontaneous symmetry breaking in bilayer graphene 
causes changes in the low energy band structure, namely the 
position of the pockets and even the number of pockets can al- 
ter. The symmetry breaking induced by the electron-electron 
interaction in bilayer graphene can adequately be described 
by the Hamiltonian suggested by Lemonik et al. [7]. Very re- 
cently, the same form of the Hamiltonian has been derived 
by Mucha-Kruczyiiski, Aleiner and Fal'ko for electrons in 
strained bilayer graphene [23]. They studied the band struc- 
ture topology and Landau level spectrum for strained bilayer 
graphene. The Hamiltonian depends on a complex parame- 
ter w and its change causes transition in the electronic band 
structure. However, the effect of the symmetry breaking on 
the minimal conductivity is an open question. 

In this work, we calculate the minimal conductivity as a 
function of the parameter w. We find that the change of this 
parameter can dramatically affect on the minimal conductiv- 
ity. Starting from the Kubo formula we develop a general and 
simple method to find the minimal conductivity for a wide 
class of Hamiltonians. Using this general approach we derive 
an analytical expression for the conductivity tensor in bilayer 



graphene with symmetry breaking by complex parameter w. 
As a self check we performed numerical calculations based 
on the Landauer formula and the agreement is very good. Our 
analysis of the minimal conductivity presented in this work 
was inspired by recent experiment and insightful discussions 
with Novoselov |24]. 

General approach. — To calculate the conductivity of var- 
ious two-dimensional electronic systems with electron-hole 
symmetry, we consider a general model Hamiltonian: 



^(P) 



h{Y>) 
/i*(p) 



(1) 



where /i(p) = hi (p)+i/i2 (p) and * denotes the complex con- 
jugation. Specifically, the Hamiltonian for bilayer graphene 
including symmetry breaking has this form with 



h{p) = -£l 



pi 



P+ 



(2) 



valid for the valley K. Here el = mv'^ /2 is the Lifshitz en- 
ergy (where m and vs are given in Ref. yj), the dimensionless 
momentums are p± — Px i: ipy (in units of mv^). Finally, 
w is generally a complex parameter (independent of k) and 
can be originated from the electron-electron interaction [7J 
and/or from the applied strain 12311 and/or from the slide of 
thelaver HIIsll . 

To find the conductivity for clean and bulk systems we start 
from the general Kubo formula presented in Ref. 21 (derived 
from the form given by Ryu et al. 1 18]). The minimal dc con- 
ductivity (at zero frequency uj, at zero temperature and at zero 
Fermi energy) is given by 
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and p = tik, l,m = x,y, and Ud is the degeneracy (for bilayer 
graphene Ud — nsUy, where Ug = 2 is the spin degeneracy 
and rii, = 2 is the valley degeneracy corresponding to the val- 
ley K and K'). The parameter ?/ can be interpreted physically 
as a finite inverse lifetime induced by impurities. 

The eigenvalues of the Hamiltonian ([T]) are E(k) = 
±^/i*(k)/i(k). Since H^ = E^h (here h is the 2 by 2 
unit matrix), the operator inverse in Eq. (l3cb can be written 
as [if + i?^(k)] "^ = [t]'^ + E'^{k)] "^ /2. Hence, the main 
contribution of the integrand in (l3b] i comes from the vicinity 
of the zeros of the energy eigenvalues. The real solutions of 
-E(k) = are denoted by k*^*', where s = 1, ... no and no 
is the number of zeros (for bilayer graphene no = 2, 3, 4 de- 
pending on the parameter w). Expanding the matrix Tim in 



(l3c] i around one of the zeros k^*^ ie, for fixed s the denomi- 
nator becomes a polynomial of q, where q = k — k'*^). If the 
energy dispersion -E(k) at the point k'^^ forms a Dirac cone 
then the leading term of the expansion is quadratic in q (oth- 
erwise one needs to expand E'^ (k) beyond the quadratic terms 
in q), and can be written as E'^{q) = J^i ^'^ij Qi Qj^ where 
the matrix M is positive semi-definite. For the Hamiltonian 
([U one finds 
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^'' dki dkj dki dkj ' 
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where the derivations are evaluated at k = k''*^. Moreover, 



we find that Tr 



dHjk) ag(k) 
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= 2Mim and Thnisi^rf) = 
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Now substituting this expression into ( [3b] i and re-scaling 
the wave number q as q = riq' the integral Ii„i becomes 
independent of -q, and reads 
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Note that for the limit yy — > the main contribution in the in- 
tegral Iim comes from the vicinity of each Dirac points; there- 
fore, the integral over q' can be extended to infinity. 

If the determinant of the matrix M at k^^^ is zero then the 
integral is divergent. For finite determinant of M the integral 
in (|5]i, ie, the contribution to the integral dJbl i over the vicinity 
of k^''' can be performed analytically and it becomes I^'^' — 



As) 



lin 



^ M^^/VdetM(^. Here M,^^/ denotes the matrix defined 
in Eq. Q evaluated at k = k^^\ 

Then using Eq. ( l3al ) the minimal conductivity is the sum of 
the contributions from each Dirac cones at k*^*) : 
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This universal procedure to find the minimal conductivity for 
clean systems can be applied for all Hamiltonians given by 
Eq. ([T]l except for cases when the determinant of the matrix 
M is zero at any of the Dirac points. In this case, we cannot 
derive an universal and rj independent expression for the min- 
imal conductivity and one needs to evaluate Eq. ^ numeri- 
cally for finite ?/. However, if the value of rj in Eq. (O is less 
than the energy scale for which the quadratic expansion of the 
dispersion relation is valid then our expression (|6]l gives the 
same result as that obtained numerically from Eq. ([3]). In such 
cases the minimal conductivity is independent of the value of 
rj, ie, the microscopic details of the systems. 

Applications. — We now apply our method for calculat- 
ing the minimal conductivity in different models of bilayer 
graphene. 

(i) The simplest case is when the trigonal warping is absent, 
ie, the linear term in momentum is missing in Eq. (|2|i. Then 
we find that (t™"^{w) — 8 0-0^2 and independent of w and rj. 



(ii) Now we take into account the effect of the trigonal 
warping using the Hamiltonian ([T]!-©. Then, it is easy to 
calculate the matrix M defined in Eq. (|4]i and we have 



M(k) 






ZiKii 
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(7) 



and dct M = (k^ ^ 1) , where the wave number k is in units 
of mv^/h and the matrix M is in units of {hei^/mv^) . Thus, 
the singular points (when det M ~ 0) are located on a unit 
circle in the k-plane. In this case the parameter w lies on 
the triangular like curve ABC on the complex w-plane shown 
in Fig [T^ (the same is plotted in Fig. 3a in Ref. l23h . At 
these points two Dirac cones collide and annihilate resulting 
in a topological phase transition in the energy band disper- 
sion. There aie four/two Dirac points in the momentum plane 
for the parameter w lying inside/outside of the triangular like 
curve ABC, respectively. 

In what follows it is useful to parameterize the complex 
parameter w in region I and II shown in Fig[T^ in the following 
way 



-i4a 



2cos(2/3) 



,i2a 



(8) 



where a G [0,7r/3] and < /3 < min{3a,7r — 3a}] in re- 
gion I, while in region II /? is a pure imaginary number such 
that i/3 < 0. The triangular like curve ABC in Fig [T^ in 
this parametrization reads as w/\ (a) = e^*"*" + 2 e'^" with 
< a < TT. The parameter w lying outside regions I and II 
can be folded back to either of these regions by a symmetry 
operation belonging to the group C^v . Note that if no sym- 
metry breaking is present then w = which corresponds to 
a = P = 7r/6. 

Using the parametrization dH) the minimal conductivity can 
be analytically obtained from Eqs. (|6]l and (|7]i for arbitrary 
complex parameter w (for the location of the Dirac points see 
the Supplementary Information). In region I we have 



-min 



h + 
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8(7o '^ 2 sin /3 (sin 3 a — sin /3) 

cos^ a — sin /3 sin a cos a (sin /3 — cos a) 



cos a (sin /3 — cos a) sinQ!(sin/3 + sina) 
while in region II it is given by 
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cos^ /3 — cos^ 3 a 
cos^ 2 a sin 2a cos 2a 
sin 2a cos 2a sin^ 2 a 



(9b) 



The eigenvalues of the conductivity tensor: erf 2/(8170 ) = 

1 + 2/ ( 1 ± \J\ — 4sin/3 (sin 3a — sin/3) j in region I and 

af V(8c^o) = 1, while cr|V(8CTo) = 1 + 1/ (cos 2^ - cos 6a) 
in region II. The minimal conductivity for w lying outside re- 
gions I and II can be obtained by cr™*" [w) = R^ ^ ■ cr™*" [R ■ 
w) ■ R, where i? is a symmetry operation of the symmetry 
group Csi, (reflection or 120° rotation) which transforms w 



into region I or II. Therefore, the minimal conductivity has 
a C-sv symmetry in the w-plane. The minimal conductiv- 
ity tensor cr™" is a symmetric matrix but generally for com- 
plex w the non-diagonal elements can be different from zero 
and even can be negative. Note that for w — (without 
symmetry breaking) we recover our earlier results, namely 
CT™(u; = 0) = CT™"(u; = 0) = 24cro and the off-diagonal 
elements are zero |,21.] . 

We numerically calculated the components of the minimal 
conductivity tensor from (|9]l (or from Eqs. (|6]l and (|7]l) for 
complex parameter w and plotted in Fig.[lJ), c and d, respec- 
tively. One can see from the figures that the points w/\ along 
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FIG. 1. (Color online) a) The complex w-plane with the triangle like 
curve ABC wa (red solid line) and regions I and II. The dashed and 
solid lines in regions I and II correspond to constant parameters a 
and /3, respectively. The components of the conductivity tensor (in 
units of (To) in the complex w-plane: b) axxiw), c) ayy{'w) and d) 
axy{w) = Gyx{w). Here w is in units of sl- 



the curve ABC, ie, where the electronic topological transition 
occurs the conductivity changes dramatically. We compared 
our theoretical prediction (|9]l with that obtained numerically 
from ([3]) for finite 77 and found that the numerical results be- 
comes better with decreasing 77 (for details see Fig.[ST|in the 
Supplementary Information). The analytical result (|9]l starts 
to deviate from the numerical calculations at some values of 
parameter w for which the Dirac points are close enough that 
the energy corresponding to the saddle point between them 
becomes the same order of magnitude as the inverse lifetime 

To confirm our analytical predictions ^ we calculated the 
conductance using the two-terminal Landauer formula (for de- 
tails see the Supplementary Information). In these calcula- 
tions the orientation n — (cos 0, sin 0) of the graphene sample 
with respect to the contacts is fixed (see Fig.[S2]in the Supple- 
mentary Information). From the conductivity tensor (|9]l the 
conductivity along the direction n in the two-terminal mea- 



surement is given by aKubo{w) — n • (t{'w) ■ n [18, 26]. Then 
this resuhs should be compared with that obtained numeri- 
cally from the Landauer formula. Our results are shown in 
Fig. 121 The complex parameter w is taken along two different 
lines given by the direction oj in Fig[T^ such that w — wqc^'^, 
where wq varies. Similarly, two electrode directions n are 
taken. As can be seen from the Fig.|2]l the agreement between 
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FIG. 2. (Color online) The conductivity crKubo('!i'o) (black solid line) 
and by numerical evaluation of the Landauer formula (red dashed 
line) as functions of the parameter wo, where w — uioe'", and a) 
t^ = 27°, e = -42° and h) lj = 57° , 9 = 42°. The units are the 
same as in Fig.[T] The dotted vertical lines indicate the values of the 
parameter uiq where the number of Dirac points changes. 

the two approach is very good for most of the values of wa 
(and similarly good agreements were found for other values 
of w and directions n not presented here). Note that the elec- 
tronic topological transition occurs at two values of wq, where 
the number of Dirac points changes from four to two (see the 
dotted lines in Fig.|2|. Near the singular points the deviation 
between the numerical and analytical results arises from the 
fact that here the Dirac cones come closer to each other than 
the momentum space resolution of the numerical method dic- 
tated by the finite size of the sample in the Landauer approach. 
Thus, not too close to the singular points w/s our central re- 
sult ^ well approximates the elements of the conductivity 
tensor for bulk and clean bilayer graphene which is indepen- 
dent of the inverse lifetime 77 and depends on the topology of 
the Fermi surface at the neutrality point. 

We now comment the experimental feasibility of measuring 
the minimal conductivity. On the one hand, in recent exper- 
iments 1,24] the bilayer graphene samples are clean enough 
that it is ballistic up to lengths 2 — 4 fira making sense 
to apply the Landauer approaches and to compare with our 
analytical result (|9]) obtained from the Kubo formula. Fur- 
thermore, the temperature experimentally can be as low as 
T w 0.25 K w 0.02 meV, ie, its effect can be neglected in first 
approximation. On the other hand, the conductivity depends 
very sensitively on the orientation of the sample and the values 
of the parameter w. Thus, these two unknown factors seem to 
be crucial to reproduce the measured minimal conductivity. 
However, with the experimental control of the direction of the 
electrodes the measurement of the minimal conductivity may 
provide a good tool for determining the complex parameter w 
and exploring its origin in the electronic topological transition. 

In conclusion, using the Kubo formula we calculated ana- 
lytically the minimal conductivity in bilayer graphene taking 



into account the electronic topological transition. Our results 
are confirmed by numerical calculations of the conductivity 
using the Landauer approach. We hope that our analysis may 
provide a better insight into the origin of the reconstruction of 
the electronic spectrum observed in recent experiment |l24J]. 
Note that our general approach for calculating the minimal 
conductivity can be applied to various other electronic sys- 
tems. 
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DIRAC POINTS 

In the parametrization of the parameter w given by Eq. dHJ 
in the main text the zeros k'*' in region I are given by 



k(i) 



k(2) 



COS 2 a + 2 cos /3 cos a 
sin 2 a — 2 cos /? sin a 

cos 2 a — 2 cos /3 cos a 
sin 2 a + 2 cos /3 sin a 



(SI a) 



(Sib) 



where T denotes the transpose operation. For u = no 
symmetry breaking occurs in bilayer graphene, and the ze- 
ros at k*^^^ = (0, 0) and the remaining three ones are lo- 
cated at the corner of a regular triangle as well known. For 
the singular point u — —1 the zeros k^'^) and k*^^^ collide at 
k = (1,0), while for m = 3 the zeros k(i),k(2) and k^^) 
collide at k= (-1,0). 

Evaluating the matrix M at the zeros k^^^ (s = 1, . . . , no) 
and using Eq. (|6]l in the main text we find (after a simple alge- 
bra) that the universal minimal conductivity in clean bilayer 
graphene as a function of the parameter u is given by 
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Cxy (u) — for all values of u, ie, the con- 



and aa:yiu) 

ductivity tensor er is a diagonal matrix. Note that for u < — 1 

or u > 3, ie, where only two Dirac points exist ^"""(u) is 

constant. 



k(3) 



k(4) 



■ COS 2 a + 2 sin f3 sin a 
sin 2 a + 2 sin /3 cos a 

cos 2 a — 2 sin /3 sin a 
sin 2 a — 2 sin /3 cos a 



(Sic) 



(Sid) 



while in region II there are only two Dirac points and they are 
the same as Eqs. dSlab and ( ISlbb but with purely imaginary 

/3. 



REAL PARAMETER w 

For real parameter w, ie, for w = u, where w € K we can 
give analytical expressions for the components of the minimal 
conductivity tensor directly in terms of the real parameter u 
(without using the parameters a and f3). In this case, the zeros 
k'**) as functions of u are given by 



k(i) = (-1, -V3~uf, if u < 3, 



k(2) = (-l,^ 



if li < 3, 



k^^) = (1 - VuTT, of, if M > -1, 

k('*) = (i + /;7TT,o)^, ifM>-i, 



(S2a) 
(S2b) 
(S2c) 
(S2d) 



FINITE INVERSE LIFETIME 

For real parameter u our analytical result ( IS3I ) is compered 
with the numerical results obtained from Q in the main text 
and plotted in Fig.[ST] We found that the agreement between 
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FIG. S 1 . (Color online) The diagonal components of the minimal 
conductivity tensor as functions of the real parameter u obtained 
from ([S3} (solid lines) and from numerical evaluation of Eq. l[3j in 
the main text (dashed lines) corresponding to 77 = 0.09. The units 
are the same as in Fig.[T]in the main text. For — 1 < it < 3 there are 
four Dirac cones, while for u < — lorii>3 only two Dirac cones 
exist. 




FIG. S2. Schematic of a wide bilayer strip of length L, contacted by 
two electrodes. Solid lines in the lattice represent the upper layer, 
while the dashed lines the lower layer. The electrodes are at an angle 
6 with respect to the symmetry axis of the bilayer lattice denoted 
by n. The wave vector kx (ky) is perpendicular (parallel) to the 
electrodes. The carrier concentration around the neutrality point can 
be tuned by a separate gate electrode not shown in the figure. 



the theoretical predictions ( IS3l l and the numerical results be- 
comes better with decreasing rj. The reason for the deviation 
between the analytical ( |S3t and numerical results is discussed 
in the main text. 



LANDAUER'S APPROACH FOR CALCULATING THE 
CONDUCTIVITY 

Landauer's approach, in a similar way as in Refs. |l)- 
[SO, gives the conductance in terms of transmission proba- 
bilities for an electron to be transmitted by a sample. In 
our two-terminal calculations, we employ a model of infi- 
nite wide bilayer strip of length L with carrier concentration 
near the charge neutrality point as the scattering region, and 
two highly doped regions as the electrodes oriented at angle 
9 with respect to the symmetry axis of the bilayer lattice (see 
Fig. |S2t . According to recent experiments ^ in very clean 
samples we can consider the bilayer strip to be ballistic up to 
lengths 2 — 4 /im. In our numerical calculations we chose 
i = 2.5 /im. The connected electrodes are modeled by semi- 
infinite bilayer sheets under negative bias voltage Vb much 
higher than any other characteristic energy scale in the sys- 
tem. Since properties of the transmission of the strip saturate 
for Vq < ^1 eV, in our numerical calculations we optionally 
applied Vb — —10 eV. To describe low energy excitations 
around the charge neutrality point in the vicinity of the Bril- 
louin zone corners K (^ = 1) and K' (i^ = —1), according to 
Refs. |5,,6J, we consider the Hamiltonian: 



i/4(k) 



where fc 
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= e=F**(fcj^ 


± ik||) (here fcj_/fc|| is 


the com 


po- 



trodes), v — 3/2 djo/h is the Dirac velocity in the monolayer, 
^3 = 3/2 dj^/h, and d is the bond length between the carbon 
atoms, and finally w is the same parameter than in the 2x2 
Hamiltonian given by Eqs. ([T]i and Q in the main text. Here 
we take 70 ^ 3 cV, 71 ^ 0.4 cV and 73 ~ 0.3 eV for the 
intra- and interlayer hopping parameters, respectively 151]. In 
bulk systems the low energy band structure is well approxi- 
mated by a two by two effective Hamiltonian when the high 
energy dimer states are decimated |J7|]. However, in this case 
this decimation cannot be used directly since a potential step 
generates further terms in the effective Hamiltonian. 

The band structure around K and K' are identical up to a 
rotation by vr in momentum space. Since the conductance is 
invariant under this transformation, it is sufficient to perform 
the calculations only at K point, since the other valley would 
lead to the same contribution to the conductance. 

The transmission probabilities is obtained by solving the 
scattering problem. Electronic states are specified by their en- 
ergy e and the transverse momentum fc|| , which are conserved 
in the scattering process. For a given e and fc|| there are four 
solutions of longitudinal wave vector k± which satisfies the 
characteristic equation det [iJ4(fcj^, fc||) — /4e] = 0, where 

I4 is the 4 X 4 identity matrix. Solutions fcj"' (n e {1, . . . , 4}) 
can be obtained if one rearranges the characteristic polyno- 
mial into the form: 



det [Hi{kj_,k\\) - he] = Yl «»»(£, fcji)fcj 



(S5) 



m=0 



The complex roots k\^ of the polynomial ( IS5I ) are equal 
within numerical precision to the eigenvalues of the compan- 
ion matrix of this polynomial ISJ. Electronic states in the 
scattering region (0 < 2: < L) are denoted by '^^c (x, y) — 
<j)(n)gi(fc^ K+fcyy)^ where <l>s"' satisfy relation for all possible 
n: 



F4(fcr,A:||)cl>L"^=e1>L"^. 



(S6) 



The scattering state in the bilayer strip is then a linear com- 
bination of these four electronic states. Similarly, in the elec- 
trodes (x < or x > L), one can obtain the longitudinal 

in) 

wave numbers k\^ \^^^ and the corresponding electronic states 



,(") 



n:d(a^,y) = *& 



- ^i- 



lead' 



^(fc'r.Ldaj+feiiw) 



using Eqs. (ISSl ) and ( IS6b 
with substitution e — > £ — Vq. Generally, inside each electrode 
two states are propagating or decaying to the left (labeled by 
L) and other two to the right (labeled by R). If we assume 
two incident states in the left electrode, than the ith (i = 1, 2) 
scattering state takes form: 



vI/W(x,j/) 



nent of the wave vector perpendicular/parallel to the elec- 



n=l 


'■«*Ltd" 


'(- 


y), 


x < 0, 


n=l 








0< x <L 
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where r*j and i*j are the reflection and transmission ampli- 
tudes, which have to be determined (together with coefficients 
An) by imposing the continuity condition of the wave func- 
tions at the boundaries x — and x = L. Then the conduc- 
tivity, including valley and spin degeneracies, reads yl]: 



two Dirac cones align in the direction kj^, the resulting con- 
ductivity will show quantum interference of electronic states 
of the two cones. 



^Landauer 
CTO 



2L / dfci 



EEK 



(S8) 



where co ~ e^ /{ttK). In the numerical calculations we find 
that the conductivity a is practically independent of L in the 
ballistic regime with lengths L = 2 — 4 jiva. 

Finally, we note that (i) the conductivity calculated from 
Eq. ( IS8I ) becomes sensitive to peculiar dynamics of the ne- 
matic phase transition when the length of the bilayer strip {L) 
is much larger than the characteristic length of the band struc- 
ture I ^ 7r/Afc, where Afc is the splitting of the Dirac cones. 
For w = the characteristic length of the band structure is 
I ^ 50 nm ||31], which is much smaller than the length of the 
bilayer strip L — 2.5 /im used in our calculations, (ii) when 
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